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The Chandrasekhar mass limit sets the scale for the late evolutionary stages of massive stars, 
including the formation of neutron stars in core collapse supernovae. Because its value depends on 
the gravitational constant G, the masses of these neutron stars retain a record of past values of 
G. Using Bayesian statistical techniques, I show that measurements of the masses of young and 
old neutron stars in pulsar binaries limit G/G = (—0.6 ± 2.0) x 10 _12 yr _1 (68% confidence) or 
G/G = (-0.6 ± 4.2) x 10" 12 yr _1 (95% confidence). 
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a. Introduction Whether or not the fundamental "constants" of nature vary with time has been a question of 
\& ' considerable interest since Dirac suggested that the gravitational force may be weakening with the expansion of the 
universe Although general relativity predicts G is identically zero, a variable G is expected in theories such as the 
Brans-Dicke scalar-tensor theory and its extensions , and has recently received renewed attention in the context 
of extended inflationary cosmology Q . 

The most direct experimental limits on G/G come from monitoring the separations of orbiting bodies, since from 
Kepler's laws it is easily shown that G/G — —a/ a, where a is the orbital semimajor axis. Early studies based on 
ancient occultation and eclipse observations had relatively low precision but strong limits are now available 

from direct measurements. For convenience, I define C-12 = (G/G)/10~ 12 yr~ 1 . Lunar laser ranging experiments yield 
C-12 = ± 11 |(J, while radar ranging to Viking gives C-12 = 2 ± 4 |ll[] or —2 ± 10 p2]| , depending on assumptions 
about solar system mass uncertainties and correlations between model parameters. Similarly, observations of the 
pulsar-white-dwarf binary PSRB1855+09 yield C-12 = — 9 ± 18 (u|. (Observations of the double neutron-star 
binary PSRB1913+16 give C-12 = 4± 5 fli)] , but this limit is greatly weakened when the G-driven variation in the 
r*** ' gravitational self-energy of the stars is considered fll5|| .) 

Indirect evidence about past values of G can be obtained from comparison of big-bang nucleosynthesis models with 
the observed 4 He abundance ]l6| . A recent reanalysis argues that 0.7Gq < Gbbn < l-4Go 17|, corresponding to 
IC-12I 0.9 for a power law variation of G, or IC-12I ^ 40 for a linear variation of G. Limits are also obtained 
from considerations of the long term stability of clusters of galaxies and globular clusters [jl8| (IC-12I ^ 40 — 60), or 
from evolution of the Sun or other stars, since the luminosity of a star L oc G 7 jl9],^0]]. Unfortunately, the Earth 
preserves only a crude memory of the early luminosity of the Sun, so paleontological evidence gives only a weak limit, 

h : ic-12] < 100 

Fortunately, a much more precise record of early stellar evolution can be found in the galactic population of neutron 
stars, whose masses are set at their time of formation by the balance between the Fermi degeneracy pressure of a cold 
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electron gas and the gravitational force, through the Chandrasekhar limit |2 



M <* ~ \g^ n ) ' (1) 

where mjv is the mass of the neutron. Because M c h sets the mass scale in the late stages of stellar evolution [M , we 
expect the average neutron star mass fi ~ M c h, which implies G/G = —2/i/3[i. In this Letter, I show how observations 
of neutron star masses and ages can be used to set tight limits on /i and hence G. 

b. The Masses and Ages of Neutron Stars There are now five double neutron star binaries known. In each case, 
five Keplerian parameters can be very precisely measured by pulse timing techniques p4[ : the binary period Pb, the 
the projection of the orbital semimajor axis on the line of sight x = aisini, the eccentricity e, and the time and 
longitude of periastron, To and o^o- These parameters are related to the pulsar and companion masses, mi and m2, 
through the mass function 

_ (m 2 sku) 3 4tt 2 x 3 
; ~ (m!+m 2 ) 2 ~ ~GPj- ( ' 

In each case, the relativistic advance of the angle of periastron, d>, has also been measured, which yields an estimate 
of the total system mass m t = m x + m 2 . For three systems (PSRs B1534+12, B1913+16, and B2127+11C) the 
measurement of the combined effects of the transverse Doppler shift and the gravitational redshift allow the individual 
determination of the pulsar and companion masses. In Table 1, I collect the measurements of /, mi, m.2, and mt, 
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and their uncertainties (the uncertainty in / can be neglected). To a good approximation, the measurements arc 
independent and normally distributed, so if rhi is the true value of one of mi, 7712, or m, t , then the probability that 
we will measure mi is 

An upper limit to the age of a pulsar can be determined from the rate at which the pulsar period P is increasing 
because of the loss of rotational energy to radiation, t c = P/2P. The progenitors of the double neutron star systems 
were binary systems consisting of two massive stars that underwent successive supernova explosions. Because the 
lifetime of a neutron star progenitor is only about 10 7 yrs, the age difference between the pulsar and companion can 
be neglected. In the case of PSR J1518+4904, for which t c is at least 16 Gyr, I assume instead an upper age limit of 
10 Gyr, about the age of the Galactic disk. 

PSR B2127+11C is a special case, because it is in the globular cluster NGC 7078 (M15), where interactions can 
result in new companions being exchanged into binaries. Indeed, in the standard scenario, neutron star formation is 
completed within the first ~ 10 7 yrs of the cluster lifetime, so the small characteristic age of the pulsar is understood 
as the time since the pulsar was (last) spun-up by a companion. In this case, the relevant age for both neutron stars 
is the cluster age, which can be found by comparing stellar structure models with observations of the color-magnitude 
distribution of cluster members. The most recent calculations, using the latest nuclear equation of state and opacity 



data, find a cluster age of 12-13 Gyr |25|], somewhat younger than previous estimates |26[|27|]. 

c. Variability of the Average Neutron Star Mass In Figure 1, I display the available information on the ages and 
masses of neutron stars in double neutron star binaries. There is certainly no evidence that the average mass has 
changed by more than a few tenths of a solar mass in the last 12 Gyr. I introduce a model in which the average mass 
varies as /i = /io — /it, where t is the neutron star age, with the goal of estimating the posterier density P(fi\{xi}, {ti}), 
where {x{\ and {ti} are the observations of neutron star masses and ages. 

The underlying distribution of neutron star masses is unknown, but the tight clustering of masses of young stars, 
for which G is unimportant, suggest that a normal distribution with variance s about fj, is reasonable. Hence the 
probability that a neutron star of age t will have mass m is 

P(m\t, Mo , A, a) = -±= e-C^-w-^Va.. (4) 



I consider a normal distribution because the present data set is too small to justify a more complex model; other 
distributions (e.g., uniform between an upper and lower bound |p8fl ) are possible, but the final results are fairly 
insensitive to the form of Eq. (|]). 

Assuming a uniform prior density for (1, I use Bayes' Theorem to write P(fi\{xi}, {ti}) oc P({xi}\(i, {ti}), where the 
proportionality constant is set by the normalization condition J P(fi\{xi},{U})dfi = 1. Then 

P({Xi}\il,{ti}) 

P({xi}\fi, fi , s, {U}) 7r(// ) 7r(s) d/j, ds, (5) 



where I take the prior density n(no) to be uniform for positive [Iq, and tt(s) as uniform in logs |29| ]. 

To evaluate Eq. (||), note that the independence of the measurements of the five binary systems allows the factor- 
ization 

P({xi}\fi,fi ,s,{ti}) = JJP(xj|/i,/x ,Mi) = II J P(xi\fj,fiQ,s,ii)P(ii\ti)dii, (6) 

i i 

where now Xl = {m 2 ,m t } for PSRs B1534+12, B1913+16, and B2127+11C, and x t = {f,m t } for PSRs J1518+4904 
and B2303+46. I take P(U\ti) uniform for < ti and zero otherwise, except for PSR J1518+4904, where I make 
the further (conservative) assumption that P(ti\ti) = for ti > 10 Gyr, and for PSR B2127+11C, for which I take 
P(ti\ti) = S(ti — 12Gyr). We can further factor 

P(m 2 ,m t \fi,fio,s,ti)= / P(m 2 \m 2 )P(m t \mt)P(m2\li,Vo,s,i)P(m t -m c \i2,^o,s,i)dmtdm2, (7) 



and a similar, more complex expression for P(f, m,t\fi, (J,, s, U), using Eq. (|^) and P(cosi) = 1/2. Using Eqs. (ffl and 
(§), we can then (numerically) evaluate P(fi\{xi},{ti}). I find /i = -1.2 ± 4.0(±8.5) x 10" 3 A/ Q Gyr _1 at the 68% 
(95%) confidence level, corresponding to G/G = -0.6 ± 2.0(±4.2) x 10~ 12 . 
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d. Discussion The measurement C-12 = —0.6 ± 2.0 is a factor five tighter than earlier limits. It is important to 
understand how sensitive this result is to our model assumptions, the most critical of which concern the neutron star 
ages. Independent evidence that pulsars with small spin-down rates are old comes from studies of pulsars in binaries 
with cooling white dwarf companions. For example, an upper limit on the optical luminosity of the companion of 
PSR B1855+09 yields a minimum system age of 4Gyr, comparable to the pulsar timing age t c = 5Gyr J30| , ^3[ . 

Some post- formation mass transfer is required to recycle the pulsars to their observed periods, however a very 
small transfer (UOAMq) is sufficient in all cases (and the physics of spin-up is presumed independent of time). Of 
more concern is the age of PSR B2127+11C, which is based on the standard model of neutron star formation: core 
collapse in massive stars. It has been proposed that some neutron stars in globular clusters may result from accretion 
induced collapse (AIC) of O-Ne-Mg white dwarfs 31 32]], though theoretical and observational uncertainties remain 
p2fl . In the most unfavorable case, both PSR B2127+11C and its companion were (separately) formed by AIC before 
exchanging into the current binary. The pulsar may then be young, with ti < t c , while the companion could have any 
age less than the cluster age. Using these assumptions, I find the weaker limit C-12 = 2.3 ± 5.0. 

Mass determinations of pulsars such as PSR B1855+09, for which white dwarf cooling ages are available, are clearly 
of great interest. In fact, observations of the Shapiro time delay in this system yield mi = 1.501q '^Mq Jll|, but the 
uncertainties are still too large for this measurement to contribute significantly to our estimate of C-12. 

Time variability of the Chandrasekhar limit has other potentially observable implications, most notably for Type la 
supernovae, which are widely believed to be the thermonuclear disruptions of Chandrasekhar mass white dwarfs |j3"| , 
and which are remarkable standard candles, with an intrinsic dispersion of only 0.12 magnitudes when multicolor light 
curve shape corrections are done f34j| . The optical emission is due to decay of 56 Ni. In a naive model, 56 Ni production 
will be roughly proportional to mass, and hence to G~ 3 ^ 2 . If G is proportional to the expansion parameter R a , then 
G/G = <jHq. An ambitious program on existing telescopes could measure the average luminosity of supernovae at 
z = 1 to ~ 0.05 magnitudes [§5[, corresponding (for Hq = 60kms _1 Mpc _1 ) to the magnitude change produced by 
C-12 = 2.7. 
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FIG. 1. Masses of the neutron stars in the five binaries of Table 1. Circles indicate individual stars; squares are the average 
mass in cases where the individual masses cannot be determined. Ages shown are upper limits, except for PSR B2127+11C, 
as described in text. Ages of the components of a binary are offset slightly for clarity; uncertainties in the mass of a pulsar and 
its companion are not independent. The variation in the average neutron star mass corresponding to f-12 = ±10 is shown. 
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